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1 Introduction 



In recent years lots of progress has been made in QCD lattice calculations. One important progress in 
the light quark sector concerns the values of the quark masses that can now be reached. These are very 
close to the physical ones making a controlled, i.e. trustable, chiral extrapolation of the lattice results 
to the physical points possible. A very powerful model-independent framework to perform this ex- 
trapolation is Chiral Perturbation Theory (ChPT), the Effective Field Theory of QCD at low energies. 
Indeed it allows to calculate low-energy QCD processes in terms of the light pseudoscalar mesons 
masses. Hadron properties are presently actively studied on the lattice and chiral extrapolations to 
their physical values are performed, see for example HI El El . 

At the same time lots of effort is put into testing the Standard Model (SM). In order to do so one 
has to have very precise determinations of the QCD quantities which generally enter the different 
processes under consideration. Two very interesting quantities in this respect are the strangeness 
changing scalar f and vector / + form factors which are measured in fQ 3 decays [0. Indeed a 
measurement of the K e3 inclusive decay rate leads to the extraction of the product of the vector 
form factor at zero momentum transfer /+(0) and of the CKM matrix element \ V US \. Consequently 
the knowledge of /+(0) allows to extract this matrix element and thus to test the unitarity relation 
between the elements of the first row of this matrix. Another test comes from the values of this form 
factor at the Callan-Treiman point and at its soft-kaon analog @. Indeed, at these particular 
points the scalar form factor has a well-known value as dictated by SU(Nf) x SU(Nf) low energy 
theorems, with Nf = 2 and Nf = 3, respectively. Combining this information with experimental 
results from semi-leptonic decays one can determine the values of the scalar form factor at these two 
points in the SM. Thus a departure from these values would be a sign for physics beyond the SM 
such as right-handed quark couplings to the [|5l 13 or charged Higgs effects, see for example the 
discussion in [4J(and references therein) and (S). However, in order to have a reliable and accurate 
test of the SM one should know very precisely the corrections to the Callan-Treiman theorem and its 
soft kaon analog which are only exact in the soft meson limit. They are usually calculated in ChPT 
BH. In Ref. [6J the one-loop result from Ref. [9] was used and an estimation of the higher order effects 
was done since at next-to-leading order some low-energy constants (LECs) contribute which are not 
very precisely known at present. Experimentally there has recently been interest in trying to obtain 
the value of the scalar form factor at the Callan-Treiman point. The three collaborations NA48 IfTOl , 
KLOE ||TT| and KTeV |[T2l have reanalysed their data so as to extract this value using in their analysis 
a dispersive representation of the form factors proposed in Refs. ||5l[6]. With the current experimental 
precision the data from the last two collaborations show a good/marginal agreement with the SM 
while NA48 has a 4. 5 a deviation. 

The scalar form factor has been studied on the lattice. Some parameterization of its momentum- 
dependence plus the knowledge of the one-loop ChPT result at zero momentum transfer is used to 
extract /+(0). Here we will fit the lattice data from Ref. [13] for the scalar form factor using a ChPT 
calculation at two loop order 031. Furthermore, we will not only consider the scalar form factor 
but at the same time we will fit the ratio of the kaon to the pion decay constants Fx/F n lfT5l [T6l 
since, as we will see, similar LECs enter the two quantities. This will allow us to determine some 
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LECs at two-loop order (0(p 6 )) and thus not only obtain /+(0) and determine \V US \ but also the 
value of the scalar form factor at the Callan-Treiman point and at its soft-kaon analog. Of course 
one should keep in mind that we are dealing here with SU(3) quantities which involve the strange 
quark mass. The question is whether one should consider the strange quark as light compared to the 
QCD scale A ~ 200 MeV or should it be treated as heavy. Related to that is the question whether 
standard SU (3) ChPT which assumes that the quark condensate is large, is a well converging series, 
the relevant expansion parameter being in that case (mx/A x ) 2 ~ 0.4 2 . Also ss sea quark pairs may 
play a significant role in chiral dynamics leading to different patterns of chiral symmetry breaking in 
Nf = 2 and Nf = 3 chiral limits lfT71[T8l . For example, lattice QCD seems to indicate a problem in 
the extrapolation of Fk/F^ to its physical value when using SU(3) ChPT to one loop order |fT6ll while 
a fit within "Kaon ChPT" lfT9]| where the kaon is treated as a heavy particle leads to good agreement. 
The mass dependence of the scalar form factor has been studied within this scheme in Ref. [20J. We 
will use here standard ChPT to two loops and we will study the convergence of the chiral expansion. 
We will also discuss the leading order 0(p 4 ) LEC LI which is related to the Okubo-Zweig-Iizuka 
(OZI) rule violation. 

In section [2l we discuss briefly the scalar form factor at two loops in ChPT. We present the lattice 
calculations in section [3] and discuss our fits and results in section HI We conclude in section |5] 



2 ChPT to two loops 

The strangeness changing form factors are defined from the K — > n matrix element of the vector 
current = s^^u 

(^(p^)\s^ t ,u\K(p K )) = (pk+Pk)^ f+(t) + (PK-Pn)^ /-(*)> (2.1) 

where t = q 2 = (px ~ Pn) 2 - The vector form factor /+(£) represents the P-wave projection of 
the crossed channel matrix element (0|s7 At w|A'7r) whereas the S-wave projection is described by the 
scalar form factor defined as 

Ht) = f + (t) + 2 * a /_(t) . (2.2) 

At zero momentum one has 

/o(0) = /+(0) . (2.3) 

These form factors were calculated to two loops in ChPT in Ref. lfT4ll . These authors introduced the 
quantity 

kit) = f + (t) + , * ^ (f-(t) + 1 - it) = fo(t) + ^ 2 * f 1 - ■ ( 2 - 4 ) 



m 2 K — mlV v " 7 ' F n J JlJV " / ' m l K — m\\ F w 

The two-loop expressions of the two decay constants F K and F w can be found in Ref. lfT5l . They 
involve two Lj, L r 4 and L\, at 0(p 4 ) and the four 0(p 6 ) Ci, C[ 4 , C[ 5 , C[ 6 and C[ 7 . Assuming that the 
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LEC L 4 is small, which is in principle the case in the standard scenario of ChPT, one can expand, as 
usually done, the denominator in the ratio of the two decay constants so that its contribution to order 
p 4 cancels and one is left with the contribution from L5 and two combinations of three as detailed 
below: 



F K /F n = l + ±( m 2 K -ml)Ll + — 2 n,l(C\, \ 2(" n ) (2.5) 

7T ^0 



2m\m\ ( -(C[ 4 + C[ 5 ) + \{C[ b + 2C[ 7 ) ) + 2m A K {C u + C 1Sj 



+ 5 



where F is the pion decay constant in the chiral limit. We will come back to the discussion of this 
equation in Section 0] 5 contains the loops and the contributions of the Lj at 0(p 6 ). Interestingly 
the dependence on these LECs is exactly the same in fo(t), see Ref. |fl4|. Thus the main advantage 
in considering f is that this quantity has no dependence on the L\ at order p 4 , only via order p 6 
contributions and furthermore, it only depends on the two 0(p 6 ) LECs C[ 2 and Cg 4 . Its explicit 
dependence on those is given by 

/o(*) = l-4(^2 + C , 34)K-^) 2 + 8^(2C[ 2 + ^ 4 )K + m2) 
^0 ^0 

-^t 2 C[ 2 + A(t) + A(0), (2.6) 

where we used the notations of Ref. [14J. As before, the quantities A(t) and A(0) have contributions 
from loops and from the LECs L\ at 0(p 6 ) and can in principle be calculated to order p 6 accuracy 
with the knowledge of the L[ to order p 4 . Parameterizations of these quantities in the physical region 
of Kis decays can be found in Ref. [fl~4|l . 

Eq. (|2.4I) is in fact inspired by the Callan-Treiman theorem [|2TT| which predicts the value of f (t) 
at the so called Callan-Treiman point, t = Ak-k = rn 2 K — m\ in the SU(2) x SU(2) chiral limit. One 
has 

/ (A^) = ^ + A CT , (2.7) 

where A C t is a correction of O It has been estimated within ChPT at next-to-leading order 

(NLO) in the isospin limit flU with the result 

A^° = (-3.5 ±8.0) • 1(T 3 , (2.8) 

where the error is a conservative estimate assuming some typical corrections of 0(m u> d) and 0(m s ). 
From Eq. (12.61) one can calculate the contribution from the 0(p 6 ) LECs to Act- It reads 



16 

A CT \ Ci = — 2 {2C{ 2 + Cl,)m\{m\ - m\) . (2.9) 
^0 



0(p 6 ) calculations [|22l using some estimates for the LECs C[ 2 and Cg 4 give results consistent with 
Eq. (12.81) . Strong isospin breaking as well as electromagnetic effects have also been evaluated Il22ll23l . 
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Another interesting quantity is the soft-kaon analog of the Callan-Treiman theorem (|24| 

f (A Kn ) = ^ + A CT , (2.10) 

with A KlT = —A K7T . A one loop calculation of the SU(3) correction A C t m the isospin limit BUl 
gives Act = 0.03. This is larger than its soft-pion analog A C t, see Eq. (12.81) . by a factor m 2 K /m 2 ^, 
however, rather small for a first order SU(3) x SU(3) breaking effect, which is expected to be of the 
order of about 25%. 

The value of V u d, the first element of the CKM matrix is very accurately known from superallowed 
+ — > + nuclear /3-decays [|25l 

\V ud \ = 0.97425 ± 0.00022 . (2.11) 

Combining this value with the experimental value of the branching ratio I" * k V2(i) /Fn l2(y) and as- 
suming the standard couplings of quarks to the W-boson allows to determine the ratio of the decay 
constants Fk/Ft,. Using instead the inclusive decay rate r# ie3 , . 11261 . one obtains the value of 
the vector form factor at zero momentum transfer /+(0). From these information and Eqs. (|2.7I) . 
(|2.10l) . one can deduce the value of the normalized form factor at the Callan Treiman point C = 
f (A KlT ) / f + (0) and at A Kn . For the explicit formulae and more details see for example Ref. I1261 . 
One has the following updated values in the SM 

f+(0)\sM = 0.959 ±0.005, (2.12) 

Fk/F v \ sm = 1-192 ±0.006, 

In C\ SM = 0.2169 ±0.0034 + A OT //+(0) , 

f {A K7r )/f + {0)\ SM = 0.8302 ± 0.0074 + A CT // + (0) . 

Deviations from these SM predictions would thus be a sign of new physics. For example at NLO 
within the minimal not-quite decoupling electroweak low-energy effective theory (LEET) Il27ll . in 
the light quark sector one has two combinations of parameters of spurionic origin describing the 
couplings of quarks to the VT-boson to be determined from experiment [|28l l6l. While the knowledge 
of the scalar form factor at the CT point measures one combination, its knowledge at Ak-k measures 
the other one. A precise determination of A C t and A C t would thus help to settle the issue of the 
presence of right-handed couplings of quarks to the W -boson. 

In order to have a very precise determination of /+(0) as well as A C t an d A cr , one needs to have 
a very precise determination of all the LECs L\ and CI which enter Eqs. (|2.5U2.6I) . 

• The Li have been determined in Ref. [|29l from a fit to the masses and to i^ 4 -decay data from the 
E865 experiment, assuming that L\ and Lg are 1/N C suppressed and using F K /F W = 1.22 0. 
Matching the dispersive results for the subthreshold expansion parameters of tiK scattering 

3 In this fit some of the Cj are taken from resonance saturation, the others are set to zero, see Ref. 1031 
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Fit 10 D291 


nK Roy Steiner (H 


Prelim. Fit All(*) 051 


Lattice l[T6ll 






set a 




set b 








0.432 


1.05 ±0.12 


0.99 ±0.13 


— 


10 3 L r 2 




0.735 


1.32 ±0.03 


0.60 ±0.21 


— 


10 3 L r 3 




-2.35 


-4.53 ±0.14 


-3.08 ±0.47 


— 


10 3 L r 4 







0.53 ±0.39 


0.70 ±0.66 


0.33(0.13) 


10 3 L r 5 




0.97 


3.19 ±2.40 


0.56 ±0.11 


0.93(0.073) 


10 3 L r 6 









0.14 ±0.70 




10 3 L r 7 




-0.31 




-0.21 ±0.15 




10 3 L T 8 




0.6 




0.38 ±0.17 




10 3 (2L£ 
10 3 (2L£ 


-LI) 








0.032 (0.062) 


-LI) 








0.050(0.043) 



Table 1: 0(p 4 ) LECs at a scale \i = 0.77 GeV. 



with their chiral expansion at order p 4 ll30l leads to somewhat different results, especially L\ is 
suggestive of a significant violation of the OZI rule in the scalar sector, see Table [T] This is in 
agreement with a determination of some of the LECs in an analysis of J ftp decays into vector 
mesons and two pseudoscalars Oil . 

• In Ref. Il32l it was shown that it was possible to reproduce the values of the Li in terms of 
properties of the light meson resonances (masses and coupling constants). The idea of using 
resonance saturation also for the 0(p 6 ) LECs was thus taken up and the Cj are presently mostly 
estimated in that framework Il33ll34l l35l . There are, however, a few problems. First the scale 
at which they are obtained is not known. It is usually assumed to be given by the lightest 
scalar nonet that survives in the large iV c limit, Ms = 1.48 GeV. The value at another scale, 
typically the p mass scale, is obtained using renormalization group equations. Furthermore, a 
test of the naturalness of the Cj 11361 shows that some of them are in fact not dominated by the 
resonance contributions. Also the LECs we are interested in have important contributions from 
the scalar sector where one knows that the OZI rule is strongly violated and where the presence 
of the wide scalar a and k mesons makes the calculation in terms of tree level diagrams from 
a resonance Lagrangian not really appropriate. Considering more specifically C\ 2 and C^, 
several calculations have been performed based on the study of the scalar form factor with 
AS = [37] or AS = 1 [f38l|3ll. In the literature these two LECs, Eq. $ZB, lie in the range 
-lO- 3 GeV" 2 to a few 10" 4 GeV~ 2 . The four other 0(p & ) LECs (C[ 4 , C{ 5 C[ 6 and C[ 7 ), 
Eq. (12.51) . needed in our study are not very well known. In Ref. |[40l where the Ci have been 
recently determined within a quark model, one finds C[ 5 = C[ 6 = 0, C[ 7 = 0.01 • 10~ 3 GeV -2 
and C[ 4 = —0.83 • 10~ 3 GeV~ 2 which is smaller than what is found in resonance saturation 
CL = -4.3 ■ 10~ 3 GeV" 2 . 



With the progress of lattice QCD it becomes also possible to extract the LECs from a chiral extrapo- 
lation of the lattice data. Already some of the 0(p 4 ) ones have been obtained mostly within SU(2) 
(k). Relations between the SU(2) and the SU(3) LECs allows to determine the L { from the k [|4Tll42l 
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(for similar relations between the C{ see Ref. 114310 . Results from the RBC/UKQCD collaboration are 
shown in Tabled] As can be seen from this table most of the 0(p 4 ) LECs are still not well enough 
determined for a very precise test of the SM. A global fit of all the low-energy constants of Chiral 
Perturbation theory at next-to-next-to-leading order currently performed P4l l45l will hopefully help 
to settle the values of these LECs much more precisely. Some preliminary results P31 which differ 
from fit 10 by using some more recent data, by letting L\ and U 6 free and by adding some constraints 
from nK scattering show better agreement with the analysis of Ref. [30] as the comparison between 
the second and third column of Table Q] shows. 



3 Lattice 



Following the pioneering work of Ref. [46 ] different collaborations have extracted the vector form 
factor at zero momentum transfer either with N f = 2 07] HI g9] |50l or N f = 2 + 1 [Q21 flavours. 
The idea is to first evaluate the scalar form factor f (t) at the momentum transfer t max = {m K — m 1 ,)' 1 . 
This can be very efficiently done calculating a double ratio of three-point correlation functions [46] . 
Then a phenomenologically motivated interpolation is performed up to zero momentum transfer^ and 
the Ademollo-Gatto theorem is used to obtain a rather precise value for /+(0). Let us consider the 
chiral expansion of /+(0) 

/+(0) = l + / 2 + /4 + --- , (3-1) 

where f n = C((mx A /(47rF 7r )) ri ) and the first term is equal to one due to gauge invariance. The 
Ademollo-Gatto theorem BTI states that the deviation from unity of /+(0) is predicted to be second 
order in SU(3) symmetry breaking, i.e. of order (m s — m) 2 , where m s and rh are the strange and 
the average of the u, d quark masses, respectively u so that the 0(p 2 ) term f 2 in the chiral expansion 
of /+(0) is free of any LECs. The different collaborations generally take this term from a one-loop 
ChPT calculation [H 

h = -0.0227 , (3.2) 

obtained for pion, kaon and eta masses taken at their physical values and in the isospin limit and 
determine the difference 

A/ = / + (0)-l-/ 2 . (3.3) 

This difference contains of course all terms starting at the order 0(p 6 ). Also used is the partially 
quenched expression derived in Ref. Il53~1 . An expression for f 2 using NLO SU(2) ChPT can be 
found in Ref. |]20l . The first determination of A/ in a quark model framework gave Af = —0.016(8) 

EH. 

The RBC/UKQCD collaboration for example JO UM simulates with N f = 2 + 1 flavors of 
dynamical domain wall quarks. In order to determine /+(0), they performed a simultaneous fit to 

4 A new technique has been developed in Ref. [52] which will allow to directly simulate at t = on the lattice. 
5 Note, however, that despite this theorem the light quark mass difference m u ^ mg can modify /+(0) to first order. 
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both the t and quark mass dependences using the ansatz 



fo(t, m v , m K ) = — — — . (3.4) 

1 - t/(M + Mi{m z K + m%)Y 

This formula motivated by the Ademollo-Gatto theorem has four fit parameters A , Ai, M , Mi and 
f 2 is the NLO term, Eq. (13.31) . They have also used a second order Taylor expansion as parameteriza- 
tion of the t-dependence of the form factor @. They obtain 

/+(0) = 0.9644 ± 0.0033 stat ± 0.0037 syst . (3.5) 

The same collaboration has also extracted the ratio Fk/ F n lfl~6ll 

F K /F n = 1.205 ± 0.018 stat ± 0.062 syst . (3.6) 

A summary of other lattice results can be found in Refs. Il2~6l 1551 . In the unquenched simulations 
they fall in the range between 1.189 and 1.218 for the central value of Fk/F^ and between 0.956 and 
0.968 for the one of /+(0). While the errors on the former are very small, they are larger on the latter. 
All these numbers should be compared to the Standard Model values, Eq. (12.121) . 



4 Chiral Extrapolation 

We now turn to the central point of the paper, namely the chiral extrapolation of the lattice data 
on Fk/Ft, and /+(0) based on the two-loop ChPT calculations |[i"4||. We use the results from the 
RBC/UKQCD collaboration since this is the only collaboration which has calculated both these quan- 
tities with Nf = 2 + 1 flavors. We take the data performed on the 24 3 x 64 volume with an inverse 
lattice spacing of a -1 = 1.73(3) GeV and a simulated strange quark mass, am s = 0.04 close to 
its physical value. We do not correct for finite volume effects (FV) or lattice artefacts (LA). They 
have been estimated for Fx / F n lfi~6l where the error bars they quote for these effects are roughly 
equal (FV) or even larger (LA) than the statistical ones. We only included the statistical errors in 
our fits. Also we did not include the correlations between Fk and F n since they are not available. 
Lattice results have been obtained for four values of the light quark masses which correspond to pion 
(first number in parenthesis) and kaon masses (second number) equal to (0.329, 0.575) GeV (set (I)) 
(0.416, 0.604) GeV (set (II)), (0.556, 0.663) GeV (set (III)) and (0.671, 0.719) GeV. Clearly, ChPT 
cannot be valid at too high pion and kaon masses so we completely discard the last set in our fits 
and mostly use sets (I) and (II). For each pion mass they have calculated the scalar form factor at 
five values of t going from ~ —0.4 GeV 2 to t max - Again for the fits we only use the three smallest 
absolute values of t. 

6 This parameterization and the pole one are usually assumed either in lattice calculations or in most of the experimental 
analyses. One should note, however, that the pole parameterization has no real physical motivation in the case of the scalar 
form factor. Also it has been shown lfT2l l6l that in order to get a very precise parameterization of the scalar form factor in 
the physical region of decay (m| < t < t max ), an expansion up to third order had to be done. 
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Fit I 


Fit II 


Fit III 


Fit IV 


Fit V 


Fit VI 


Cl2 


5.77 ±0.56 


7.84 ± 0.58 


4.69 ±0.95 


5.74 ±0.95 


4.69 ±0.56 


4.43 ± 0.88 


C*34 


2.54 ±0.43 


-1.28 ±0.44 


3.76 ±0.95 


1.07 ±0.96 


3.76 ± 0.43 


3.50 ±0.94 


C14 


0* 


0* 


0.65 ± 1.38 


0.71 ± 1.42 


0.65* 


-0.93 ±0.67 


2C 17 


0* 


0* 


0.31 ±3.31 


1.92 ±3.36 


0.31* 


4.16 ± 1.56 


F 


89.8 ±0.1 


69.2 ±0.0 


89.8 ±0.1 


69.3 ±0.0 


89.8* 


89.8 ±0.1 


/+(o) 


0.956 


0.963 


0.956 


0.961 


0.956 


0.958 




1.20 


1.19 


1.20 


1.19 


1.20 


1.19 


InC 


0.22 


0.20 


0.22 


0.21 


0.22 


0.21 


/o(A^) 
10 3 A cr 


0.75 


0.75 


0.75 


0.76 


0.75 


0.77 


1.00 


-2.14 


0.27 


-3.65 


0.18 


-0.32 


10 2 A CT 


-9.00 


-9.86 


-8.24 


-8.18 


-8.11 


-7.03 


10 3 A 


18.08 


17.77 


18.24 


17.66 


18.18 


16.71 


x 2 


1.40/4 


0.96/4 


1.67/4 


1.29/4 


3.01/4 


4.8/7 



Table 2: Result of the fits to the lattice data. The first five quantities are the parameters of the fits. The 
star denotes an input quantity. The Cj's are in units 10~ 4 GeV 2 and F is in MeV In Fits (I) and (II) 
the lattice data on F K / F T are not included. 



A quantity O at two loops has typically the following form after renormalization of the pion decay 
constant 

0(m m m K , m v ) = U LO ± H =j— , (4.1) 

where Olo, C^nlo an d Onnlo are the contribution at leading order (LO), next- to-leading order (NLO) 
and next-to-next-to-leading order (NNLO), respectively. Here, F w is the pion decay constant calcu- 
lated at 0(p 4 ) at the value of the pion mass and of the kaon mass under consideration and F is the 
pion decay constant in the chiral SU(3) limit. When working at the physical pion and kaon masses 
one usually replaces everywhere the decay constant by its physical value, since the difference is of 
higher order. This is the procedure which has been used to determine the LECs, set (a) and (b) of 
Tabled] If this is mostly justified for set (a) where the difference between F w and F is small, this is 
more questionable for set (b) where F = 67.1 MeV but allows, of course, for a better convergence 
of the chiral series. Also going away from the physical point the difference between F n and F might 
become again too large for this procedure to be entirely satisfying. Here we will just replace F by the 
physical value of F n in the NNLO term in order to be consistent with the determination of the LECs. 
Also to be consistent with their determination we will use Eq. (|2.5I) for determining Fk/F^. Again if 
this is justified for set (a) where the convergence of this quantity is rather good as we will see below, 
this is more questionable for set (b). In the expression Eq. (14.11) . the mass of the r] enters the NLO 
and NNLO terms. In the calculation of Cnlo its NLO expression is used while in Onnlo the rj mass 
is given by the Gell-Mann-Okubo relation. 

We have performed several fits to the lattice data and determined from these fits results for /+(0), 
the slope of the scalar form-factor at zero momentum transfer Ao, Fx/F n , Act and Act- We have 
taken F as a parameter of the fit using the value of the physical pion decay constant as input. Apart 
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from Fit (VI) they are done with the two lattice data sets with the smallest pion values, sets (I) and 

(II) . All the fits are done for the three smallest absolute values of t. The values of the LECs Li are 
taken from sets (a) and (b) of Table [TJ These sets correspond to a value of m s /rh = 24. In Ref. Il44ll 
another preliminary set is given corresponding to a somewhat larger value m s /rh = 27.8 as obtained 
by MILC and HPQCD/UKQCD. It leads to an even smaller value of F = 62.7 MeV and will not be 
discussed here. The results of the fits are given in Table|2] The one from this other set are comparable 
to the one of set (b). 

• Fits (I) and (II) are three parameter fits of f (t) using sets (a) and (b) respectively. The Cj are the 
one used in the determination of the 0(p 4 ) LECs, fit 10, namely C[ 4 = C[ 5 = C[ 6 = C[ 7 = 0. 
For set (a) Fk /F n = 1.22 whereas for set (b) Fk / F n = 1.19. Slightly different values are given 
in the table for set (a) since, as explained below Eq. (|4.1I) we did not use the physical value of 
Fk in the calculation of this quantity in the NLO term but rather its NLO expression. 

• Fits (III) and (IV) are combined fits of F K /F W and f (t) using sets (a) and (b), respectively, 
as in the previous fits but now the combinations C[ 4 + C[ 5 and C{ 5 + 2C[ 7 which appear in 
Fk/ are left free. Since we need to determine F , we, in principle, need to know C[ 6 and the 
combination C[ 5 — 2C[ 6 . We will assume them equal to zero, this is consistent with the results 
in Ref. [|40l . Thus we do in fact determine C[ 4 and C[ 7 . 

• Fit (V): here we fix the combinations C[ 4 + C[ 5 and C[ 5 + 2C[ 7 from Fit (III) and we fit the 
quantity /„(*). 

• Fit (VI) is the same as Fit (III) but with the lattice data for f (t) from set (III) also included. 

As can be seen from Tabled we obtain very good fits of the lattice data. Fits (I) and (II), however, 
do not reproduce well the two lattice points for Fk/F^ from sets (I) and (II). Fits (III) and (IV) which 
correspond to two very different values of L\ are comparably good, but an order of magnitude larger 
value of C{ 7 is in fact needed in order to compensate for the larger value of L 4 in Fit (IV) compared to 
Fit (III). C[ 4 and C[ 7 are at least an order of magnitude smaller than what is expected from resonance 
saturation in the scalar sector which leads to typical values ~ 1CT 3 . One has for example [TT3TI 

^ CjfirrAn ^ 3 Q y-2 (4 2) 

where q, c m and d m are coefficients of the scalar chiral Lagrangian. M s and d m are obtained from 
the masses of the scalars 7^(1430) and a (980) and c m = 0.042 GeV and c d = 0.032 GeV The 
results of Fit (VI) do not differ much from Fit (III), only Ci 4 and Cyj are larger in absolute value and 
the slope of the scalar form factor is somewhat smaller. This fit is shown on Fig. Q] for sets (I) and 

(III) . Even though we only fit the three smallest points in absolute value, the t-dependence of set (III) 
is remarkably well reproduced by ChPT to two loops. 

Fitting f (t) leads to strong anticorrelations between C 12 and C 34 on the one hand and C 14 + Ci 5 
and C*i5 + 2C*i 7 on the other one, typically of the order of —0.8 while in Fit (V) the correlations 
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t [GeV~; 




Figure 1 : Momentum-dependence of the scalar form factor. The results of Fit (VI) are displayed 
(solid line) for set (I) (left panel) and set (III) (right panel). The convergence of the chiral expansion 
is also displayed: the dash-dotted line is the result at LO, the dashed line displays the one up to NLO. 



between C12 and C34 are reduced by a factor of two. Also a comparison of Fits (III) and (V) shows 
that the error bars on these two LECs are smaller in the latter case. Thus a rather good determination 
of the LECs Ci 2 and C 34 is obtained by fitting the function f . Their order of magnitude is the one 
expected from resonance saturation. Note that the value obtained for Ci 2 + C 34 is rather independent 
of the fits within one set, one gets ~ 8 • 10~ 4 GeV~ 2 for set (a) and ~ 6 • 10~ 4 GeV -2 for set (b). 

The results for Fk/F^, /+(0) and InC are consistent with the values obtained assuming the 
standard quark couplings to the W -boson, Eq. (12.121) . We refrain to give error bars here since one 
should have a more precise knowledge of the Li as well as lattice data at lower pion and kaon masses 
to really be able to pin down these quantities very precisely. Difference between the various sets gives 
an idea of the errors. The value of A turns out to be rather large compared to the experimental results, 
the lattice determination of Ref. Il50ll or to what is obtained from the formula obtained in a dispersive 
parameterization of the form factor [51 

2 

TfL 

A = - - w — (In C- (7(0)), G(0) = 0.0398 ±0.0044, (4.3) 

where G{t) is a dispersive integral of the phase of the form factor which is identified in the elastic 
region with the s-wave, 1 = 1/2 Kw scattering phase according to Watson's theorem. In the analysis 
[21 it was taken from [3(3 where a matching of the solution of the Roy-Steiner equations with the 
Kw — y K-k , 7T7T — > KK and Tin — > ixix scattering data available at higher energies has been per- 
formed. Note that in this analysis the LECs obtained, second column of Tabled! are more consistent 
with the values used in Fit (IV), especially a large violation of the OZI rule was found. This large 
value of Ao can be traced back to the too large value of the combination 2C\2 + C34 which enters its 
expression within ChPT, see Eq. (12.61) . It is however compensated by a small curvature X' Q leading to a 
value of the scalar form factor at the Callan-Treiman point in agreement with the SM value. Typically 
one obtains X' ~ 1 ■ 10~ 4 instead of ~ 6 ■ 10~ 4 as expected from experiments and dispersive analyses 
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Il56l l5l. Again Ci 2 has a too large positive value. Stringent constraints on slope and curvature have 
recently been obtained using the method of unitarity bounds 11571 . 



Let us study the convergence of the results. In Fig. [T] is shown f (t) as obtained in Fit (VI) at 
LO (dot dashed line), NLO (dashed line) and NNLO (full line). On the left-hand-side (LHS) set (I) 
is displayed and on the right-hand- side (RHS) set (III), in order to compare the dependence on the 
pion and the kaon masses. Clearly, as expected, the convergence of f (t) worsens as one increases 
the absolute value of t (LHS, set (I)) and as one increases and m K . At the physical pion and kaon 
masses one has from Fit (III), 

/ + (0) = 1 - 0.019 - 0.026 + ..., 
F K /F n = 1 + 0.140 + 0.061 + ..., 
fo(A Kn ) = 1 + 0.139 + 0.063 + ..., (4.4) 
A CT = - 0.0025 + 0.0028 + ..., 
Act = + 0.024 -0.106 + ..., 

and from Fit (IV) 

/ + (0) = 1-0.019-0.019 + ..., 
F K /F n = 1 + 0.113 + 0.081 + ..., 
f (A K7T ) = 1 + 0.110 + 0.081 + ..., 
Act = - 0.0033 - 0.0003 + ..., 

Act = + 0.021 - 0.103 + ..., (4.5) 

where the first, second and third terms are the 0(p 2 ), 0(p A ) and 0(p 6 ) contributions, respectively, 
and the ellipses denote terms of order p 8 and higher. Note that by definition Act and Act have no LO 
contribution. The convergence is rather good/not very good for Fx/F n and fo(A.Kn) for set (a) and 
set (b) respectively while the one for A C t is good for set (b) and not for set (a). One should however 
keep in mind that the NLO correction for this last quantity is small being an SU (2) x SU (2) one. Also 
the NNLO contribution of A C t is of the expected size of the corrections, Eq. (12. 8b ■ The convergence 
of /+(0) and Act is bad whatever the set. However the convergence looks again worth than it is 
in reality. Indeed for both quantities the contribution at NLO is smaller than naively expected. For 
/+(0) this is essentially due to the Ademollo-Gatto theorem, as we have seen in the previous section. 
In both cases the NNLO term is of the expected size. Let us look in a bit more details at the diverse 
contributions for /+(0), F K /F n and f (A K7r ). One has for Fit (III) 



/+(0) 


= 1 + 


(- 


■0.019 + 0.000) + (0 


.012-0.003 


- 0.034) 


+ 


Fk/ F n 


= 1 + 


(0 


.057 + 0.083) + (-0 


.005 + 0.045 


+ 0.021) 


+ 


/o(A^) 


= 1 + 


(0 


.055 + 0.083) + (-0 


.001 + 0.047 


+ 0.017) 


+ 



(4.6) 
and for Fit (IV) 

/+(0) = 1 + (-0.027 + 0.008) + (0.012 - 0.002 - 0.029) + ..., 
F K /F W = 1 + (0.086 + 0.027) + (-0.005 + 0.078 + 0.009) + . . . , 
fo{A Kw ) = 1 + (0.083 + 0.026) + (-0.001 + 0.063 + 0.019) + .... (4.7) 
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The first brackets give the contribution from the loops and the Li at fourth order and the second 
brackets represent the one at sixth order from the two-loops, the one-loop with one Li insertion plus 
tree graphs with two Li and the tree graphs ~ d, in order. One sees that the large contribution of 
/+(0) at NNLO is due to big corrections of the dimension six operators, as was the case for the slope 
and the curvature, see the discussion before. It could be that the corresponding LECs Cj are larger 
than they are in nature mocking up some higher order effects. The contributions from the two-loop 
and the one-loop ~ Lj topologies do converge. In the case of Fk/F^ and fo(A.Kw) it is the terms 
proportional to Li which are responsible for their not so good convergence in the case of set (b), 
explaining the difference between the two sets. Let us consider also the convergence of F K jF^ at 
larger pion and kaon masses. One has for Fit (IV) 

F K /F W = 1 + 0.043 + 0.093 + ... = 1.136 + ..., set (I) 

= 1 + 0.023 + 0.076 + ... = 1.099 + ..., set (II). (4.8) 

For comparison the lattice data are: 

F K /F n = 1.134 + 0.011, set (I) 

= 1.101 + 0.010, set (II). (4.9) 

As already stated for the scalar form factor and as expected, the convergence gets worse when in- 
creasing the values of m, and vtik- This bad convergence could be an artefact of the use of lattice 
data obtained at still too high pion and kaon masses for ChPT to really be valid. 



5 Conclusion 

We have done here a first exploratory study using a two-loop ChPT calculation to fit the lattice data. 
Certainly finite volume effects for example should be taken into account in a more refined treatment. 
However, before this can be done, a better knowledge of the Li are necessary and more lattice data 
at smaller masses are needed. This is important for checking the convergence of the 577(3) ChPT 
calculations as well as for a more precise determination of the quantities studied here. Also if the 
result of set (a) is not very sensitive to the treatment of the NNLO term, see discussion below Eq. (|4.1I) . 
this is clearly not the case for set (b) and our results here are certainly not the final ones. Indeed, if 
large values for L\ and Lg as expected from a large violation of the OZI rule were confirmed in the 
future then the use of standard ChPT as done here would not really be appropriate. A way of solving 
the problem could be for example to work within resummed ChPT [|58ll . A study along this line is in 
progress ll59~l . 
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